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ABSTRACT: 

An  analysis  for  the  transient  axial   response  of  a  gun  launched  motor 
case  with  unbonded  propel lant  is  performed.     The  present  effort  extends  a 
previous  work  by  including  higher  harmonics  in  the  representation  of  the 
breech  force  at  the  aft  end  of  the  motor  case.     Including  the  higher  har- 
monics leads  to  a  more  accurate  representation  of  the  actual    breech  force. 
The  results  of  the  analysis  indicate  that  the  dynamic  effects  on  the  in- 
ternal  axial   force  in  the  case  are  negligible  for  a  breech  force  duration 
of  0.026  seconds.     This  is  the  same  conclusion  that  was  made  in  the  pre- 
vious work  based  upon  a  simple  representation  of  the  breech  force. 
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NOTATION 

A  area  of  motor  case  cross  section 

A,  Fourier  cosine  coefficients  of  breech  force 

B.  Fourier  sine  coefficients  of  breech  force 
c  speed  of  wave  propagation 

e  axial   strain 

E  Young's  modulus  of  elasticity 

f  axial   force 

K  constant 

I  length 

m  mass  of  motor  case 

P  breech  force 

t  time 

u  deformation 

u  displacement 

u  rigid  body  displacement 

U  transformation  variable 

v  transformation  variable 

x,  y  coordinate  axes 

X,  Y  coordinate  axes 

a  time  duration  of  breech  force 

Y  constant 

p  mass  density  of  the  motor  case  material 

a  axial   stress 

to  constant 


n 


INTRODUCTION 

This  report  examines  the  transient  axial   response  of  a  gun-launched 
rocket  motor,  with  unbonded  propellant,  during  the  launch  phase.     It  is 
part  of  a  continuing  effort  to  develop  structural   design  procedures  for 
gun-launched  rockets.     In  a  previous  report  (Ref.   1),  an  analysis  for  the 
time-dependent  axial   force  distribution  in  the  motor  case  was  performed 
in  which  the  force  at  the  breech  end  of  the  case  was  approximated  by  a 
one-half  sine  wave  over  the  duration  of  the  loading.     That  analysis  estab- 
lished that  for  the  one-half  sine  wave  input,  a  time  duration  of   .012 
seconds  leads  to  negligible  dynamic  effects,   i.e.,  for  a  one-half  sine 
wave  of   .012  seconds  duration,   the  internal   axial   force  at  any  location 
along  the  case  appears  as  a  duplicate  of  the  breech  force,  but  of  smaller 
magnitude.     It  was  shown,  however,   that  decreasing  the  time  of  the  force 
duration  to   .0012  seconds  leads  to  significant  dynamic  effects.     This  re- 
sult indicated  that  higher  harmonics  in  the  breech  force  could  cause  im- 
portant dynamic  effects.     Since  the  actual   breech  force  is  not  a  simple 
half  sine  wave,  but  in  fact  is  a  complicated  function  of  time,  an  investi- 
gation was  made  to  consider  the  problem  further. 

To  account  for  the  higher  harmonics  in  the  breech-force,   the  analysis 
decomposes  the  actual   force  into  a  Fourier  series  with  27  terms.     This 
force  is  applied  to  the  case,  and  the  internal   axial   force  is  computed  at 
several   locations  along  the  case  over  the  duration  of  the  loading.     This 
analysis  then  establishes  the  contribution  of  the  higher  harmonics  in  the 
breech  force  to  the  dynamic  stresses  in  the  bar.     The  analysis  is  restricted 
to  the  free-free  motor  case. 


ANALYSIS 


In  the  figure  below,  AB  is  the  motor  case,  treated  as  a  one- 
dimensional    bar,  under  the  applied  breech  force  P(t).     The  actual 
displacement  of  any  point  of  the  bar  with  respect  to  the  X,  Y, 
fixed  coordinate  system  is  denoted  u.     This  motion  consists  of  two 
parts,  a  rigid  body  motion  u  and  a  strain-causing  deformation  u. 


*  Y 


Then 


u  =  u  +  u 


->X 


(1) 


If  m  is  the  total  mass  of  the  bar,  then  the  rigid  body 
motion  is 


dt 


d  u     P(t) 
2     m 


The  differential  equation  of  motion  for  any  point  of  the  bar  is 

2       2 
3  u    2  3  u  _  n 

— 2  -  c  — 2  -  0 

3t^      3X 
where  c  is  the  speed  of  wave  propagation,  and  x  is  a  coordinate 

axis  parallel  to  x  moving  with  velocity  u.  The  constant  c  is 
given  by 


c  = 


where  E  is  Young's  modulus  of  elasticity  of  the  bar  material  and 

p  its  specific  gravity.     Combining  equations   (1)  and   (3)  gives 

2  2 

3        /"  \  2  3   u   _  n 

— P   (u  +  u)   -  c     — p-  =  0 

3t  3X 


(2) 


(3) 


(4) 


A 


Substituting  for  ..2  from  (2)  yields, 

a  L 

4  +    PM.C2    A    =    0> 

3t2  m  3X2 


2 


With  the  notation   (    )       denoting  3  \   '    ,  and   (   )..   denoting       [  '    , 

the  governing  differential   equation  is 

u++  -  c2u       =  -  ^^-  t  >  0,  /rx 

tt  xx  m  '  (5) 

0  <  x  <  I 

The  initial   conditions  and  boundary  conditions     for  the  free-free 
system  are 

I.C.'s:     u(x,o)   =  o  u.(x,o)   =  o  (6) 

B.C.'s:     ux(o,t)   =  $1    „xCfc,t)   =  o  (7) 

The  solution  of  equations   (5)  through   (7)  gives  the  deformation  of 
any  point  for  any  time  greater  than  zero.     The  stress  is  then  simply 
obtained  from  Hooke's  Law,  a  =  Eu   . 

The  solution  method  proceeds  by  transformation  of  the  given  system  to 
a  system  with  homogeneous  boundary  conditions.     This  is  achieved  by 
considering  a  solution  in  the  form 

u(x,t)   =  v(x,t)  +  U(x,t).  (8) 

Boundary  condition  equations   (7)   become 


vx(o,t)   =  ^-  -  Ux(o,t) 


vxU,t)  =  -uxU,t) 


(9) 


These  equations  become  homogeneous  when  U(x,t)  is  taken  as 

U(x,t)  =  (x  -  ^  )  £M  (10) 


With  equations   (8)  and   (10),  governing  equations   (5)   through   (7) 
become 

vtt  "  c2vxx  "  f<x^  (11) 

I.C.'s:     v(x,o)  =  -U(x,o)       vt(x,o)   =  -Ut(x,o)         (12) 
B.C.'s:     vx(o,t)   =    o  vt(£,t)   =  o  (13) 

where 

f(x.t)   =  -Uu  +  c2Uxx  -  ffl 

2         p  (14) 

-  (y    -    2<_    \       I"  2P_ 

^x       11  '     AE     "     PIA 
A  Fourier  series  solution  which  satisfies  boundary  conditions   (13)   is 
assumed  for  v(x,   t) ,   i.e., 

oo 

v(x.t)  -      ^vn(t)  cos  "-f  (15) 

n=o 

Accordingly,  equations   (11)  and   (12)  become 

dZ\(V  2  ,     * 

2 +     wnvn(t)   =  fn(t)         n  =  ]  '   2'    ■"  (16) 

dt 

I.C.'s:     vn(o)   =  o 

•iL    x2 


dVn   (Q,       2    /       \*~  ^!     dPW  CQS  njx  dx  (17) 

dt"  (0)   "  I   J  EA  dt  t 


where 


and 


fn(t)   =     |    /       f(x,t)  cos  ^dx,  (18) 


nirC 
%  =  ~T  (19) 


To  this  point  the  analysis  is  identical    to  that  presented  in  reference 

1.      In  that  analysis,   P(t)  was  approximated  as  a  half  sine  wave  P  sin   nt/a; 

here,   the  breech  force  is  approximated  by  the  Fourier  expansion, 

N 
PCt)   -  A0  +  £      |    Ak  cos  &S  +  ^  sin  2M   1  (20) 

k=l 

For  the  linear  system,  the  response  due  to  equation  (20)  is  obtained  by 

superposition  of  individual  terms. 

(i)  for  the  uniform  force  A  : 

o 

2P 
From  equation  (14),  f(x,t)  =  — jjr   .  Then  equation  (18)  gives 

f  =  o  for  all  n.  The  initial  conditions  are  v  =  o,  and 

v  =  o  from  equation  (17),  and  hence  v  (t)  =  o.  Then  the 
n  n 

stress  is 

a-Eux-EUx-(l  -£)£>  (2D 

This  is  the  stress  due  to  the  "rigid  body"  acceleration  for 
constant  force, 
(ii)  for  the  cosine  terms,  A^  cos  — —  : 

From  equation  (14)  we  obtain, 

l*-£)       (Zk«)2        2c2 


f(x,t)=  \-nr--  \«J  -  tfe   -Akcos  a 


Then  equations    (17)  and   (18)  give 


d 


2kfrt 


dt  vn<0)  =  °  (22) 

,  ,i>  8k  t        .  2kitt 

fn(t)  =  "  -T7Z    Ak  cos  —  (23) 

a    n   tA 


With  equations   (22)  and   (23),   the  solution  of  equation  (16) 

Vn(t)     7^2~"77  '  cos  V  )  <24) 

where 


K" = "  xfe~  (25) 

n  a  EA 


and 

2kTT 


Yn  =       ^h~     '  (26) 


Then  the  stress  associated  with  A,  cos  ^^  terms  is 

K      a 


a  -  Eux  -  E(vx  ♦  Ux) 


A 


k 1  ...  2knt 


(y'-.D 


mrx 


— —  |  cos  __  _  cos  ^t  I  •  sin  t 


+  (     |)^cos^  (27) 


The  latter  term  in  equation  (27)  is  the  contribution  associ- 
ated with  the  "rigid  body"  acceleration. 

iii)  for  sine  terms,  B.  sin  — —  : 

k  a 

Equation   (14)  gives 

x2_  2 

f(x>t)   =     J    X  "21        l2k^\  2c2     )    .  2kTit 

~Em"  \~^~j      "   EA£    J   Ak  cos     a 

From  equations   (17)  and   (18), 
,  4k£B. 

EAnSa 

a    n   EA 


The  solution  of  equation   (16)  with  equations   (28)  and 

(29)   is 

2K  kB.a2       f    y2  -  2   1  K  a2B. 

w   f  +  \   -       n     k  J    ^n c.  .    ,     n       k  .     2kTrt     ,,M 

vn(t) — j-     <    -T—      -sin  u,  t  +     2     2  sin—     (30) 

V  L  ^n  "  1    J  -  (Yn  -  1 ) 


where 

Rn  =  T^T  (31) 

The  stress  associated  with  B,    sin  =^^    is 

2 

a  =    V"       2Bk     j      2k       /  Yn  "   2   \       .         +L  1  .     2k*t).     nnx 

/  .       zrr    i     —  — o sin  aj_t  +  -o sin  fsin 


^     A  1  D  O  ol"  ^      u      '  O  ->  '  "     fain  o 

n-rrA  K  I       Z         ,      /  n  2         ,  a  <- 

n-1  L      n       V  Yn  -  1    /  Yn  ~ 

n       x\     k  «•■;„  2k^t 

+  (1    "  J]  /T  sin  T~  (32) 

Again  the  latter  term  in  equation  (32)  results  from  the  "rigid 
body"  acceleration. 

The  stress  at  any  point  is  obtained  as  the  superposition  of 
equations  (21),  (27),  and  (32). 


RESULTS  AND  DISCUSSION 


Equations  (21),  (27),  and  (32)  for  the  stress  have  been  evaluated 
on  the  NPS  -  IBM  360/67  digital  computer  for  specific  values  of  A.  and 
B.  obtained  by  a  Fourier  decomposition  of  the  given  breech  force.  This 
decomposition  is  done  by  the  IBM  subroutine  FORIT  on  the  computer.  The 
present  analysis  resolved  the  breech  force  into  26  sine  and  cosine  har- 
monics, i.e., 

(33) 


26 

I 

k=l 


P(t)=Ao+£     {     Akcos^l    +Bksin^i) 


The  magnitude  of  the  harmonic  coefficients  A.    and  B.    depend  upon  the 
nature  of  the  breech  force.     For  this  analysis,   the  force  given  in  Fig. 
1    of  Ref.   2  was  used.     Approximation  of  that  force  by  equation  (33)  leads 
to 


k 

Ak 

Bk 

k 

Ak 

Bk 

0 

14.1680 

0. 

14 

-.3164 

-.0186 

1 

-8.6082 

9.6038 

15 

-.1682 

-.0940 

2 

-3.2336 

-4.3103 

16 

-.0552 

.0146 

3 

1.3139 

-.8107 

17 

-.1202 

.0983 

4 

-.1162 

.1077 

18 

-.2053 

.0899 

5 

.1089 

-.1257 

19 

-.1821 

.0827 

6 

-.2036 

.3394 

20 

-.2162 

.0815 

7 

-.3478 

- . 0886 

21 

-.2723 

.0286 

8 

-.3001 

-.1664 

22 

-.2164 

-.0294 

9 

-.2229 

-.3071 

23 

-.1599 

.0250 

10 

.1728 

-.2293 

24 

-.2146 

.0638 

11 

.2260 

.1067 

25 

-.2589 

-.0065 

12 

-.0425 

.3243 

26 

-.2162 

.0034 

13 

-.3157 

.2287 

- 

A  plot  of  equation  (33)  with  the  above  values  of  A.  and  B.  is 
given  in  Fig.  1  for  a  .026  second  time  duration  of  breech  force. 
Figures  2  through  5  are  the  results  for  the  stress  in  the  bar  due  to 
the  breech  force  shown  in  Fig.  1  as  a  function  of  time  for  4  ratios 
of  0,  1/4,  1/2,  and  3/4,  respectively.   (The  analysis  is  for  the  5  in. 
steel  shell  of  0.19  in.  thickness).  Comparing  Fig.  i  with  Figs.  2 
through  5  reveals  that  the  internal  force  appears  as  a  duplicate  of 
the  breech  force,  with  the  stress  curves  scaled  versions  of  the  force 
input.  In  Figs.  6  through  8,  the  contributing  harmonics  are  individu- 
ally displayed.  The  elastic  dynamic  effects  are  the  small  oscillations 

about  each  harmonic  itself. 

it 

Figures  10  through  J&  are  the  results  of  the  analysis  for  a 

time  duration  of  0.0052  seconds  for  the  brech  force.     This  is  a  re- 
duction in  the  actual   duration  by  a  factor  of  five.     This  sample  prob- 
lem was  considered  in  order  to  assess  the  influence  of  the  magnitude 
of  the  force  duration.     The  elastic  dynamic  effects  can  be  clearly 
seen  in  these  figures.     Note  that  the  magnitude  of  stress  has  been 
significantly  changed  by  the  dynamic  effects. 

CONCLUSIONS 

The  results  of  this  analysis   (Figs.   1-5)  show  that  dynamic 
effects  are  negligible  for  the  breech  force  of  0.026  seconds  time 
duration  of  Figure  1,  Reference  2.     Hence,  when  the  propellant  is 
not  bonded  to  the  motor  case,   the  case  responds  to  the  breech  force 
as  if  it  were  applied  in  a  static  sense,  i.e.,  the  maximum  internal 
axial    load  in  the  case  can  be  computed  on  the  basis  of  the  mass 


distribution  and  the  maximum  rigid  body  acceleration  of  the  rocket. 
This  is  the  same  conclusion  arrived  at  in  Ref.  1  based  on  a  simple 
one-half  sine  wave  pulse;  hence,  accounting  for  the  more  complicated 
actual  breech  pressure  pulse  confirmed  the  original  conclusion. 

Further  investigation  (Figs.  10-14)  verified  the  dependence 
of  the  dynamic  effects  on  the  presence  of  significant  higher  harmonic 
terms  in  the  breech  force.  Such  harmonics  may  be  the  result  of  shor- 
tening the  time  duration  of  loading,  or  by  the  presence  of  a  local 
rapid  perturbation  of  force  on  an  otherwise  acceptable  smooth  breech 
force  input. 
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